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ABSTRACT 

We present photometry and spectropolarimetry of the pre-main sequence star 
HD 106506. A photometric rotational period of '-^1.416 ± 0.133 days has been de- 
rived using observations at Mount Kent Observatory (MKO). Spectropolarimetric 
data taken at the 3.9-m Anglo-Austrahan Telescope (AAT) were used to derive spot 
occupancy and magnetic maps of the star through the technique of Zeeman Doppler 
imaging (ZDI). The resulting brightness maps indicate that HD 106506 displays photo- 
spheric spots at all latitudes including a predominant polar spot. Azimuthal and radial 
magnetic images of this star have been derived, and a significant azimuthal magnetic 
field is indicated, in line with other active young stars. A solar-like differential rota- 
tion law was incorporated into the imaging process. Using Stokes / information the 
equatorial rotation rate, ilcq, was found to be 4.54 ±0.01 radd~^, with a photospheric 
shear 557 of 0.2ltQ Q3 radd~^. This equates to an equatorial rotation period of ^1.39 
± 0.01 days, with the equatorial region lapping the poles every ~ 30^3 days. Using 
the magnetic features, the equatorial rotation rate, 57cq, was found to be 4.51 ± 0.01 
radd~^, with a photospheric shear 50. of 0.24 ± 0.03 radd~^. This differential rotation 
is approximately 4 times that observed on the Sun. 

Key words: Line: profiles - stars: activity - stars: Zeeman Doppler imaging - stars: 
solar-type Stars - stars: individual: HD 106506 - stars: magnetic fields - stars: spots. 



1 INTRODUCTION 

Young solar-type F, G, and K stars provide proxies for the 
Sun's early evolution and insight into its intense pre-main 
sequence and zero-age main sequence magnetic activity. In 
this regard, observations of starspots and surface magnetic 
fields are important, as they provide the empirical basis for 
understanding the star's magnetic dynamo and its internal 
structure. Differential rotation is one of the primary drivers 
of the dynamo process that produces magnetic fields within 
the Sun. In the solar convection zone, differential rotation is 
constant down to the base of this zone in the region known 
as the tachochne. Strong shear forces are formed at the in- 
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terface between the solid-body rotation of the radiative zone 
and the differentially rotating convective layer. This interac- 
tion converts the poloidal magnetic field to a toroidal field. 
This effect is known as the "fi-effect". However, on young 
solar-type stars, the operation of the magnetic dynamo pro- 
cess may be occurring throughout the co nvective zone itself. 
This is k nown as a dist ributed dynamo (|Brandenburg et al.l 
1989; M oss et al.lll995l '). 

iBarnes et all l|2005l ) have investigated the link between 
differential rotation and effective temperature, and hence 
convective turnover time. They were able to fit an empirical 
power law to recent observations of a small sample of stars. 
As the temperature increases, so does the rotational shear. 
A small number of late F-/early G-type stars have been in- 
vestigated using Doppler imaging (DI) to determine their 
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Theoretical evolutionary isochrones, based on the models of Siess et al. 2000 
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Table 1. The Journal of photometric observations of HD 106506 
showing the date and filters used. The E5- Region (Graham, 1982) 
was observed on the 11th April, 2007. 
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Figure 1. The evolutionar y status of HD 106 506. This is based 
on the theoretical models of lSiess et al. I tood ) 



differ ential rotation. HD 171488 (V88 9 Her) (iMarsden et all 
l2006l ). HP 307938 (R58 in IC 2602) l|Marsden et al.ll2005l ). 
HR 1817 (lMenge]| |2005l ) and LQ Lup (RX J1508.6-4423) 
l|Donati et al.ll2000l ) have had their differential rotation de- 
termined. 

Zeeman Doppler imaging (ZDI) (jPonati et al.l l2003al ) 
enables the study of stellar magnetic fields using spectropo- 
larimetry and offers a window into the underlying magnetic 
dynamo. ZDI requires high signal-to-noise data as the po- 
larisation signa ture is typically le ss than 0.1% of the total 
light intensity (jDonati et al.|[l997l ). ZDI only measures the 
large-scale magnetic fields as small-scale magnetic fields can- 
not be recovered as the positive and negative magnetic fields 
within the resolution element are likely to be cancelled out. 
ZDI has been achieved on a small sa mple of single, G-type 



stars such as H P 171488 (V889 Her) (Marsden et al.ll2006l ) 
and HD 141943 (|Marsden etHlbOllall br 

HD 106 5 06 is a suitable target for DI and ZDI. 
Henry et al.l ( 1996 ) iden tified this star as act ive while 
Mason et all (|l998l ) and ICutispoto et al.1 (|2002l ) deduced 
this star as being single. The high resolution spectroscopic 
survey of Waitc ct al. (2005) confirmed this star as a young, 
rapidly rotating solar-type star. They measured a rotation 
rate ofiisini ~ 80 kms~^, strong Ha chromospheric activity 
(logii'Ha ~4.2), and deformation of the spectral line profiles 
indicating the presence of large starspots. 

The aim of this investigation is to map the surface of 
this star with a view to measuring the rotational shear on 
its surface and to measure and to map the magnetic field 
structure on the surface of this star. This is the third paper 
in a series that investigates the magnetic field topologies of 
young sun-like star s, with HD 141943 be ing the subject of 
two earlier papers jMarsden et al.ll2011al lbl). 



2 OBSERVATIONS AND ANALYSIS 

Spectropolarimetry of HD 106506 was obtained using the 
3.9-m Anglo- Australian Telescope (A AT) over an 11 night 
period from the 30th of March - 9th of April, 2007. Near- 
simultaneous BVR photometric data were obtained using 



UT Date 


Filters Used 


UT begin 


UT end 


2007 April 01 


B,V,R 


09:12:11 


19:03:53 


2007 April 02 


B,V,R 


11:09:03 


17:57:21 


2007 April 03 


B,V,R 


09:48:02 


18:10:43 


2007 April 04 


B,V,R 


14:33:42 


17:59:18 


2007 April 05 


B,V,R 


11:06:02 


17:14:08 


2007 April 06 


B,V,R 


09:16:22 


18:19:50 


2007 April 07 


B,V,R 


09:13:45 


18:13:12 


2007 April 08 


B,V,R 


09:06:57 


18:17:01 


2007 April 11 


B,V,R 


10:10:43 


13:26:59 



the University of Louisville's 50cm telescope at the Uni- 
versity of Southern Queensland's Mount Kent Observatory 
(MKO). MKO is located in Southern Queensland, approxi- 
mately 160km south-west of Brisbane. 



2.1 BVR Photometry at the MKO 

The University of Louisville 50cm telescope at MKO is a 
Corrected Dall-Kirkham design manufactured by Planewave 
Instruments. BVR photometric observations of HD 106506 
were taken through the standard Cousins B, V and R filters 
using an SBIG STL-630E with a Kodak KAP-6303E CCD 
chip that has an array size of 3060 x 2040, 9 fim square 
pixels givi ng a plate scale of ap proximately 0.54 arc seconds 
per pixel ||MK0 websitd 120101 . Table [T] shows the journal 
of observations while a more detailed analysis of this photo- 
metric study of HD 106506 has been included in appendix 
A. 



2.2 High Resolution Spectropolarimetric 
Observations from the AAT 

High resolution spectropolarimetric data were obtained from 
the AAT using the University College of London Echelle 
Spectr ograph (UCLES) and the Seniel Polarimeter (SEM - 
POL) (|Semellll989l : ISemel et al.1ll993l : iDonati et al.|[2003al ). 
The detector used was the deep depletion EEV2 CCD with 
2048 X 4096 13.5/im square pixels. UCLES was used with 
a 31.6 gr/mm grating covering 46 orders (orders ^ 84 to 
# 129). The central wavelength was 522.002 nm with full 
wavelength coverage from 437.71 nm to 681.56 nm. The dis- 
persion of ~0. 004958 nm at order # 129 gave a resolution 
of approximately 71000. Observations in the circular polar- 
ization signatures (Stokes V) consists of a sequence of four 
exposures. After each of the exposures, the half-wave Fres- 
nel Rhomb of the SEMPOL polarimeter is rotated between 
-|-45° and -45° so as to remove instrumental po larization sig- 
nals from the telescope and the polarimeter rDona ti et al.l 
[2003a,) . A journal of observations is shown in Table [S] 



2.2.1 Spectropolarimetric Analysis 

The initial data reduction was completed using the ESpRIT 
(Echelle Spectra Reduction: an Interactive Tool) software 
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Table 2. The Journal of spectropolarimetric observations of 
HD 106506 using the 3.9-m AAT. 



Table 3. The parameters used to produce the maximum-entropy 
image reconstruction of HD 106506, including those associated 
with surface differential rotation. 



UT Date 


UT begin 


UT end ^ 


Number 








of Sequences ^ 


2007 March 30 


09:02:56 


13:41:16 


5 


2007 March 31 


08:51:37 


13:36:05 


5 


2007 April 01 


08:51:37 


13:36:05 


5 


2007 April 02 


08:55:31 


13:38:07 


5 


2007 April 03 


09:26:01 


16:25:01 


3 


2007 April 04 


08:57:07 


12:42:40 


2 


2007 April 05 


10:13:25 


14:06:13 


4 


2007 April 06 


09:08:53 


12:07:27 


2 


2007 April 07 


08:49:55 


13:28:11 


5 


2007 April 08 


10:31:20 


11:15:43 


1 


2007 April 09 


09:35:14 


11:53:58 


2 



^ The UT end is the time when the final exposure for HD 106506 
was completed. 

^ Each sequence of 4x600 seconds exposures. After each expo- 
sure, the half-wave Fresnel Rhomb of the SEMPOL polarimeter 
is rotated between -|-45° and -45° so as to remove instrumental 
polarization signals from the telescope and the polarimeter. 



package (jPonati et al.l [l997l . l2003al ). Preliminary process- 
ing involved removing the bias and using a nightly master 
flat combining typically 20 flat field exposures. Each stellar 
spectrum was extracted and wavelength calibrated against 
a Thorium- Argon lamp. The mean pixel resolution for the 
AAT spectra was determined to be 1.689 km s~ ^ pixel" 

After using ESpRIT, Least Squares Deconvolution 
(LSD) was applied to the reduced spectra. LSD combines 
the information from many spectral lines to produce a single 
line profile, thereby provid ing an enormous mu ltiplex gain in 
the signal-to-noise (S/N) (|Donati et all 1 19971 ). The average 
S/N for a typical Stokes / LSD profile was ~1100 whereas 
the maximum S/N for the corresponding central order was 
only ~50. A G2 line mask created from the Kurucz atomic 
database and ATLAS9 atmospheric models (|Kurucj|l993l ) 
was used during the LSD process. 

In order to correct for the minor instrumental shifts 
in wavelength space, each spectrum was shifted to match 
the Stokes / LSD profile of t he telluric lin e s cont ained 
in the spectra , as w as done by iDonati et al.l (|2003al ) and 
iMarsden etHI (|2006|) Further information on LSD can be 
found in iDonati et all (|l997l '). 



3 PHOTOSPHERIC AND CHROMOSPHERIC 
FEATURES 

The brightness and magnetic im ages of HD 1 06506 were 
created using the ZDI code of iBrown et al.l (|l99ll ') and 
iDonati fc Brownl l|l997f) . It utilizes the time-series Stokes / 
and V LSD profiles as inputs for the inversion process to pro- 
duce the brightness and magnetic images. This is an ill-posed 
problem as there are an infinite number of possible solutions 
that can be found. The ZDI code uses the maximum -entropy 
minimisation techniques of lSkilling fc BrvanI (| 19841 ) to pro- 
duce an image with the mimimum amount of information 
required to produce the observed variations in the LSD pro- 
file. 



Parameter 


Value 


Photometric Period 


1.416 ± 0.133 days 


Equatorial Period (Using Stokes /) 


1.39 ± 0.01 days 


Inclination Angle 


65 ± 5° 


Projected Rotational Velocity, vsini 


79.5 ± 0.5 kms-i 


Unspotted Apparent Visual magnitude 


8.38 


Photospheric Temperature, Tp^ot 


5900 ± 50 K 


Spot Temperature, Tspot 


4000 ± 50 K 


Radial Velocity, v^ad 


13.1 ± 0.1 kms-i 


Stokes /: Qeq 


4.54 ± 0.01 radd-l 


Stokes /: dO 


0.2lto.o3 radd-i 


Stokes V: fioq 


4.51 ± 0.01 radd-l 


Stokes V: dO 


0.24 ± 0.03 radd-i 



3.1 Fundamental Parameters of HD 106506 



HD 106506 is a GIV star (|Torres et al.' '2006^1 with a me- 
dian visual magnitude of 8.5617 ([van L ceuwon 200^ . a B-V 
value of 0.605 and a V-I value of 0. 67 (based on the trans- 
formed Hipparcos and Tycho data, IESAI (Il997l')'). It has a 
trigonometric parallax of 7.96 ± 0.88 ( Ferryman et al.lll997h 
giving a distance of 126 ± 14pc or 410 ± 461yr. Using this 
V-I value, the temperature of HD 106506 is approximately 
5870 ± 45K. This was determined using the formulation of 
Bessell, Castelli & Plez (1998). The absolute magnitude is 
3.053 ± 0.244. Usi ng the bolometric corrections listed in 
iBessell et al. I 11998") and the formulation contained within 
that paper, HD 106506 is estimated to be 2.15 ± 0.26 Rq. 
The resulting luminosity of HD 106506 is 4.81 ± 1 .2 Lg . 
Figure [1] shows the theoretical tracks of ISiess et all l|2OO0i ). 
from which it may be deduced that the star is less than 10 
Myr of age and has a mass of between 1.4 and 1.6 Mq . Table 
[3]shows the initial parameters that were used in the analysis 
and the production of the Stokes / (intensity) and Stokes V 
magnetic maps. 



3.2 Image Reconstruction: Stokes / Doppler 

imaging map and surface differential rotation 

In order to produce the brightness imag e of HD 106506, 
a tw o-temperature model was used (|Collier CameronI 
Il992h . assuming one temperature was for the quiet 
photosphere, Tphot, while the other temperature is for the 
cooler spots, Tspot. Synthetic gaussian profiles were used 
to represent both the s pot temperature and the quie t 
photosphere temperature. lUnruh fc Collier CameronI (|l995l ) 
showed that there is little difference in the resulting maps 
when using synthetic gaussian profiles when compared with 
using profiles that were taken of slowly rotating standard 
stars. A number of authors have followed this formulation 
such as'Petit, Donati fc Collier Cameron] (|2002l ): |Petit et al] 
(2004) and Ma rsden et al.. (2005, 2006i ). 

The imaging code was used to establish the values of 
a number of basic parameters, including the star's inclina- 
tion angle, projected rotational velocity, t;sinj, and radial 
velocity, Vrad- This was achieved by systematically varying 
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each parameter i n orde r to minimise the y'^ value (eg see 
IJeffers fc Donatil l|2008l ): iMarsden et al.l l|2005l )'). These key 
stellar parameters were determined sequentially, first Vrad, 
theniisiui, and finally the inclination angle of the star. This 
sequence was repeated each time additional parameters were 
modified as a result of the imaging process. The full set of 
parameters that gave the minimum value are shown in 
Tableland were adopted when producing the final Doppler 
imaging map in Figure [2j a). By incorporating the photom- 
etry from MKO into the imaging process, Tphot and Tspot 
were determined by minimising the deviations between the 
measured and modelled data in both the V- and R-bands. 
The effective photospheric temperature was determined to 
be 5900 ± 50K, which is consistent with that calculated 
based on the photometri c colours of the star and using the 
bolometric corrections of iBessell et al.l (| 19981 ) within the er- 
ror bars of both determinations. The spot temperature was 
determined to be 4000 ± 50K. The photospheric-spot tem- 
perature difference of 1900K for this early G-dwarf is con- 
sistent with m e asurem en ts ma de by other authors such as 
iMarsden et al.l (l2011alf ^. 2005') and supports the relation- 
ship developed by iBerd vugina (2005). The intensity map 
with the photometry included is shown in Figure [SJ^b). 

There have been a number of diffe rent ways in which 
differential rotation have been measured. iReiners fc SchmittI 
(|2003ll have used a Fourier transform method to derive the 
parameters for rapidly rotating F-type / early G-type inac- 
tive stars. They produce a deconvolved line profile from the 
stellar spectra and then determine the ratio of the second 
and first zero of the resulting Fourier transform. This ratio is 
a measure of differential rotation on the star. However, for 
active young stars with significant asymmetry within the 
line profiles, this technique is not as effective. The differ- 
ential rotation of a star can also be estimated by tracking 
spot features at different latitudes on the surface of the star. 
A solar-like differential rotation law, as defined in equa- 
tio n [Tl has been applie d in a number of solar-type stars 
fi.e lDonati et al ] (|2003al ): [Marsden et al ] (|2005l . l2006l )). 



Surface differential rotation HD106506 



n(6i) 



(1) 



where ^{9) is the rotation rate at latitude 6 in radd^^, ficq 
is equatorial rotation rate, and dO, is the rotational shear be- 
tween the equator and the pole. In this paper the differential 
rotation of HD 106506 was measured using a technique that 
systematically adjusts the differential rotation parameters, 
fieq and dfi and determines the best fit to the data, by min- 
imising the difference between the data and the model fits 
using minimisation techniques. This technique has been 
successfully applied to a number of late F - /early-G st ars (i.e . 
Petit. Donati fc Collier CameronI (l2002l);|Pet"it et al.ll|2004 l: 



Barnes et all l|2005l ): lMarsden et alTl|2006l )). 

Initially, a fixed spot filling factor of 0.062, as deter- 
mined from the reconstructed map, was used and fleq and 
dO. were considered to be free parameters. The magnitude 
of the differential rotation was then determined by fitting a 
paraboloid to the reduced x^ values. The reduced x^ values 
for the various iloq and dO, is shown in Figure ID The as- 
sociated error ellipse, as shown in Figure |4] was determined 
by varying both the spot filling factor and the inclination 
angle. Both were changed by ± 10%. When varying the in- 
clination angle, the associated multiplications factors were 




4.51 4.52 4.53 4.54 4.55 4.56 4.57 4.58 4.59 
{rad d-') 



HD 106506. The image shows the reduced x values from the 
maximum-entropy Doppler imaging code for a fixed spot cover- 
age of 0.062 and an inclination angle of 65°. The darker regions 
correspond to lower values. The image scale projects ± 7a on 
to the axes in both Heq and dVl. The ellipse superimposed demon- 
strates the error for the value of the differential rotation and was 
determined by varying both the star's inclination and spot filling 
factors by ± 10% and determining Ccq and dTl. 



re-determined and the resulting spot filling factor as found 
to match the new angle by minimising the x^ value. 

The minimum value for f2eq was determined to be 4.54 
±0.01 radd-^and df7 = 0.211°;°^ radd" 



, The errors for the 
value of the differential rotation were determined by varying 
both the star's inclination and spot filling factors by ± 10% 
and determining fieq and dQ. respectively. By incorporating 
the differential rotation into the maximum-entropy image 
reconstruction, the resulting map is shown in Figure[2l along 
with the model fits to the LSD profiles in Figure [S] 

Figure [2] shows, on the left, the resulting maximum- 
entropy image reconstructed map generated using the spec- 
troscopic data taken at the AAT. On the right of this fig- 
ure is the reconstructed image created using both the spec- 
troscopic data and the near-simultaneous photometric data 
taken at the MKO. The photometry data has been able to 
enhance the lower-latitude features. In each map, an equato- 
rial period of 1.39 days has been used, which was established 
during the differential rotation analysis and the epoch was 
set to 2454194.93642361, which was the middle of the obser- 
vation run. 

Figure [5] shows the V and R light curves coupled with 
the maximum-entropy fits that were used to produce the en- 
hanced DI map. Figure [3] shows the model profiles overlayed 
on the measured LSD profiles. The model profiles include 
the differential rotation parameters. The dynamic spectrum 
of the residuals between the observed and modelled profiles 
is shown in Figure [S] It appears that the Doppler imaging 
process, along with incorporating the differential rotation 
parameters, has removed the majority of the large scale fea- 
tures with only a small amount of structure still apparent 
in this dynamic spectrum. 
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no photometry including photometry 

0.50 0.50 




0.00 0.00 



Figure 2. (a) The figure on the left is the maximum entropy map that was generated entirely from the AAT data, with a spot filling 
factor of 0.062. (b) The figure on the right is a similar map but including near simultaneous photometry from the MKO. The spot filling 
factor was 0.127. The value was set to 0.5 in the modelling process. The images are polar projections extending down to a latitude of 
-30°. The bold line denotes the equator and the dashed lines are +30° and +60° latitude parallels. The radial ticks indicate the phases 
at which this star was observed spectroscopically. 
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Figure 3. The maximum-entropy fits to the LSD profiles for HD 106506 with both photometry and surface differential rotation incor- 
porated into the analysis. The thick lines represent the modelled lines produced by the Doppler imaging process whereas the thin lines 
represent the actual observed LSD profiles. Each successive profile has been shifted down by 0.005 for graphical purposes. The rotational 
phases at which the observations took place are indicated to the right of each profile. 
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Reconstructed Spectra 



Figure 5. The photometric light curves for both the V filter and 
R filter (the points with error bars) coupled with the maximum- 
entropy fits (circles) for HD 106506. The bottom graph in each 
panel shows the Observed-Calculated values. 



Reconstructed Spectra 
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-50 50 

Velocity (km/s) 



100 



Figure 6. Dynamic spectra of the residuals between the 
observed-modelled LSD profiles, including differential rotation. 
The spectra is centred on the radial velocity of HD 106506. The 
three lines on the image shows the centre of the line profile (cen- 
tre line) and the rotational broadening of the line profile (outer 
lines), which is 79.5 kms~^. 



3.3 



Zeeman Doppler imaging: magnetic features 
on the surface of HD 106506 



The modelling strategy o f iDonati fc Browij (|l997D was used 
to construct the magnetic field topology on HD 106506. 
This involved utilising the spherical harmonic expansions of 
the su rface magnetic field as implemented by iDonati et al.l 
l|2006l ). The spherical harmonic expansion Imax ~ 40 was se- 



o 
a: 



-100 100 

Velocity (km/s) 

Figure 9. This figure shows a chromospheric features observed 
on HD 106506. The chromosphere is more dense at phase of ~ 
0.15. There also appears to be an area at ~ 0.75 where the chro- 
mosphere is less dense. 



lected as this was the minimum value where any further in- 
crease resulted in no further changes in the magnetic topolo- 
gies. 

Initially, a poloidal plus toroidal field was assumed, then 
the images were created by fitting the data to within the 
noise level and resulted in a mean field strength of 68.84 G. 
A similar approach to measuring the differential rotation pa- 
rameters as was used for the Stokes / (brightness) features. 
When considering these magnetic features, the equatorial 
rotation rate fleq was 4.51 ± 0.01 radd~^, with a photo- 
spheric shear 60, of 0.24 ± 0.03 radd~^. This implies that 
the equatorial rotation period was 1.38 days with a shear 
approximately 4 times that of the solar value. 

The reconstructed magnetic fields are shown in Figure 
[7] The fits to the Stokes V LSD profiles are given in Figure 

El 



3.4 Chromospheric Features: Prominences 

The Ha line has been ex tensively used as a indicato r of 
chromospheric activity (eg lThatcher fc RobinsonI l|l993l )). 

A mean Ha spectral line for HD 106506 was gener- 
ated from all the Ha profiles obtained during the observing 
run. Each individual spectrum was then subtracted from the 
mean Ha line, and dynamic spectra were then produced. 
Chromospheric features were observed at phases ~0.15 and 
~0.75. This is shown in Figure [9] 
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Figure 7. The maximum-entropy brightness and magnetic image reconstructions for HD 106506. These maps are polar projections 
extending down to -30°. The bold lines denote the equator and the dashed lines are -|-30° and -|-60° latitude parallels. The radial ticks 
indicate the phases at which this star was observed spectroscopically. The scale of the magnetic images are in Gauss. The brightness image 
(top left-hand image) has a spot filling factor of ~0.127 (or 12.7%) and is a combination of both the spectroscopic and photometric 
data. The magnetic images have a field modulus of 68.84G. Differential rotation, as measured using the Stokes V profiles, has been 
incorporated into the analysis. 



4 DISCUSSION 

4.1 Spot Features and Surface Differential 
Rotation 

Spectroscopic and photometric observations were used to 
produce tlie spot occupancy map siiown in Figure [2] Piio- 
tometry is very useful at extracting out tire lower latitude 
features (Jeffers, Barnes & Collier Cameron, 2002) but is 
not sensitive to the high-latitude features as these are al- 
ways in view, wh ereas the spectroscopy extracts t he higher 
latitude features (|Unruh fc Collier Cameron|[l997l ) but can- 
not discriminate lower latitude features at latitudes <_ 30 
degrees. Hence to get a more complete picture of the stel- 
lar image, both photometric and spectroscopic data need to 
be incorporated into the imaging process. Using the Stokes 
/ spectroscopic information coupled with the photometry, 
the spot occupancy map, as shown in Figure [2] show low- 
latitude to mid-latitude spots, and the predominant polar 



spot. This variation of spot occupancy with stellar latitude 
using spectroscopy alone and inclusive of the photometry is 
shown in Figure [TOl The fractional spottedness is defined in 
equation [2] 

F{e) = S{e)cos{6)d6/2 (2) 

where F{9) is the fractional spottedness at latitude 6, S{9) 
is the average spot occupancy at latitude 6, and d6 is the 
latitude width of each latitude ring. 

The exact reasons why rapidly rotating stars have 
polar spots are still r iot full y und e rstood. One theory 
llSchiissler et all Il996l : iBuzasil Il997l : iDeluca. Fan fc Saail 
Il997h postulates that high-latitude, but not truly polar, 
spots are a result of the increased Coriolis elfect, due to 
the rapid rotation, which may have the effect of deflect- 
ing the spot features to higher latitudes as they erupt 
thro ugh the star's convectio n zone. An alternative hypoth - 
esis (Schriiver fc Title! I2001I : iKitchatinov fc Riidigerl Il999l ) 
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Figure 8. The maximum-entropy fits to the Stokes V LSD profiles for HD 106506 with surface differential rotation, as measured using 
the Stokes V profiles, incorporated into the analysis. The thick lines represent the modelled lines produced by the Zeeman Doppler 
imaging process whereas the thin lines represent the actual observed LSD profiles. Each successive profile has been shifted down by 0.002 
for graphical purposes. The error bar on the left of each profile is plotted to ± 0.5 ct. The rotational phases, in terms of the number of 
rotational cycles, at which the observations took place are indicated to the right of each profile. 



involves the transport of flux through the process of merid- 
ional flow. It is yet to be determined whether the polar 
spots on young solar-type stars are formed at high lat- 
itudes by the strong Coriolis effect in these rapidly ro- 
tating stars, or formed at low latitudes and pushed pole- 
ward by subsurface meridional flows. Nevertheless, the issue 
can be clarified by long-term monitoring of the spot topol- 
ogy of these stars, because meridional flows should be de- 
tectable as a poleward drifts in starspo t features over time. 
IWeber. Strassmeier fc WashuettH (|2005| ) have reported ten- 
tative evidence for large poleward meridional flows on early- 
K giants, but to the best of our knowledge, such an effect 
has yet to be detected for young solar-type stars. 

These reconstructed images are consistent with those 
obtained in other Doppler images produced for young, sin- 
gle, rapidly rotating solar-type stars. The photometry en- 
abled more of the low- latitude features to b e recovered, as 
was d one for the star R58 (in IC 2602) by iMarsden et al.l 
l|2005h . 



When using the Stokes / information HD 106506 was 
found to have a photospheric shear of 0.211q qj radd~^with 
the equator lapping the poles every ~ SOlg days. When com- 
pared with the Sun's value for 5^2=0.055 radd~^, with a lap 
time of 115 days, this young pre-main sequence solar-type 
star has a shear of roughly 4 times greater than the Sun. 
This result is consis tent with the relationship determined by 
iBarnes et al. I ((2005') that differential rotation is dependent 
on spectral class (surface temperature), with M-dwarfs ro- 
tating almost as a solid body through to ear ly-G dwarfs ex- 
hibiting a st rong rotational s h ear. H owever, Marsden et al.l 
(2006) and Ijeffers fc Donatil (|2008l ) have found that the 
early G-dwarf, HD 171488 (V889 Her) has a photospheric 
shear up to seven times the solar value. So why do these 
two stars that exhibit similar surface temperature proper- 
ties have vastly different photospheric shear? One possible 
reason is that HD 106506 is a more distended star, with a 
larger convective zone when compared with HD 171488. Be- 
ing a younger star, it is probably still contracting down to 
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Fractional Spottedness versus Stellar Latitude 
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Figure 10. The fractional spottedness versus stellar latitude for 
HD106506 with and without photometry included in the imaging 
process. The fractional spottedness is based on the average spot 
occupancy at each longitude and is defined in equation [2] The 
photometry is very useful at extracting out the lower latitude 
features whereas the spectroscopy extracts the higher latitude 
features. 
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Figure 11. This shows a graph of Sil versus effective s urface tem- 
perature, as or ginallv produced bv[Barnes et al.l ll200^^;_The_data 
for HP 171488 l lMarsden et al.ll2006t ). HP 141943 l lMarsden et"aLl 
l2011al) and HP 106506 have been ad ded to this graph. H P 106506 
supports the power law developed bv lBarnes et al.l j2005 i) whereas 
HP 171488 and HP 141943 appears not to support this relation- 
ship. 



Differential Rotation using Brightness and Magnetic Features 



the main sequence. But it is still unclear why HD 171488 
should have such a large dQ and perhaps it is convective 
zone depth and the rate at which that zone is turned over 
that is more important than surface temperature (hence 
spectal type). Figure [TT] is a graph o f ^fl versus e f fective 
temperature, as orginally prod uced bvlBarnes et aL |j2005l ) 
with HD 106 506, HD 141943 (iMarsden et al.l l2011b ) and 
HD 171488 (|Marsden et all l2006l : iJeffers fc Donatil I2OO8I ) 
added to the dataset. HD 10650 6 supports the power law 
developed by iBarnes et al. I (|2005l). 



4.2 Magnetic Topology 

The reconstructed magnetic fields are shown in Figure [7] 
and the fits to the Stokes V LSD profiles are given in Fig- 
ure [S] Magnetic features in both the azimuthal and radial 
magnetic fields extend fr om low-latitu de to high-latitude 
features. As discussed by iDonati et al. I HoOSa), the radial 
component may be interpreted to be representative of the 
poloidal component while the azimuthal component may be 
considered as representive of the toroidal component of the 
large-scale dynamo field. The strong azimuthal ring around 
the star is s imilar to those observ ed on other stars such as 
HD 171488 e jeffers fc Donatill2008l ') . However, the azimuthal 
ring on HD 106506 is of different polarity to those observe d 
on other stars, such as HD 141943 (Mar sden et al.ll2011al '). 
This is unlikely to be due to instrumental polarisation as the 
data for HD 106506 were taken with the same instrumental 
setup and at the same time as HD 141943. 

The azimuthal ring of magnetic field is consis tent with 
the w reathes of magnetic fields theorised by Br own et al.l 
l|2010l) . except at a much higher latitude. Their 3- 
dimensional MHD models found that stars with increased 
rotation rate (up to 3 times that of the Sun) could produce 
self-generating wreathes of azimuthal field without the need 



Using Magnetic Features 




Using Brightness Features 
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Figure 12. This shows a graph of the error ellipses for the differ- 
ential rotation measurements for both the Stokes / (on the right) 
and the Stokes V (left). These error ellipses were generated by 
varying some of the stellar parameters, including the star's incli- 
nation (± 10°) and the Stokes / spot coverage or the Stokes V 
global magnetic field (± 10%). 



for an interface layer. One can speculate that the increased 
rotation rate of HD 106506 could be defiecting the wreathes 
of magnetic fields to much higher latitudes, as seen in Figure 
7. 

iDonati fc Landstreetl l|2009l ) found that increasing ro- 
tation rate of a star, the toroidal component of the sur- 
face magnetic field dominates over the poloidal component 
of that field. The magnetic field on HD 106506 is predomi- 
nantly toroidal with approximately 70 ± 3% of the magnetic 
energy being toroidal magnetic energy. This is similar with 
the toroidal field being the dominant field of HD 171488 
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Table 4. A comparison with other late F/early G type dwarfs that have had their surface differential rotation measured using Doppler 
imaging. 



Parameter 


IfD 106506 


HD 171488 '^'^ 


R58 (in 1C2602) 


LQ Lup ^ 


IfD 141943® 


(B-V)i 


0.605 


0.62 


0.61 


0.69 


0.7 


Photosphere 


5900 


5800 


5800 ± 100 


5750 ± 50 


5850 ± 100 


Tspot K 


4000 


4200 


3900 ± 100 


4250 


3950 


Age (Myr) 


<10 


30-50 


35 ± 5 


25 ± 10 


17 


Mass Mq 


L5 ± 0.1 


1.06 ± 0.02 


1.15 ± 0.02 


1.16 ± 0.02 


1.3 


Radius Rq 


2.15 ± 0.26 


1.15 ± 0.08 


1.18 ± 0.02 


1.22 ± 0.02 


1.6 ± 0.1 


Incl. Angle 


65 ± 5° 


60 ± 5° 


60 ± 10° 


35 ± 5° 


70 ± 10 


wsinikms"^ 


79.5 ± 0.5 


37.5 ± 0.1 


92 ± 0.5 


115 ± 1 


35.0 ± 0.5 


Stokes / Measurements 












Qcq radd"'^ 


4.54 ± 0.01 


4.786 ± 0.013 


11.165 ± 0.014 


20.28 ±0.01 


2.86 ± 0.02 


Sn radd-i 




0.402 ± 0.044 


0.0815 ± 0.026 


0.12 ±0.022 


0.24 ± 0.03 


Stokes V Measurements 










f2eq rad d~^ 


4.51 ± 0.01 


4.85 ± 0.05 






2.88 ± 0.02 












2.89 ± 0.05 


Sn radd-l 


0.24 ± 0.03 


0.47 ± 0.04 






0.36 ± 0.09 












0.45 ± 0.08 



1 Hipparcos, ^ iMarsden et al.l ||2006| ') , ^ Ijeffers h Donatil ||2008|). 
"' iMarsden et ahl | |2005^ . ' ^iDonati et al.l l l2000t). ® iMarsden et al.l l l2011al ') 
iMarsden et al.l l|2011bh 



llMarsden et al.ll2006l: iJeffers fc Donatil I2008D . AB Dor and 
LQ Hva (jDonati et aLlbOOSal ). Like HD 171488, this field 
is quite complex and much more than that operating the 
Sun. Considering the toroidal component, an octopole field 
dominates with approximately 37 ± 9% of the magnetic en- 
ergy whereas the dipole field contains only 18 ± 5% and 
a quadrople field containing approximately 16 ± 7% of the 
magnetic energy. This is in contrast to the poloidal field 
were the respective components are roughly equal between 
a dipole (10 ± 2%), quadopole (11 ± 2%) and octopole field 
(10 ± 1%). This would indicate that a much more complex 
dynamo is operating on HD 106506 and is similar to that on 
other young, sun-like stars. In addition, the toroidal field was 
predominantly axisymetric (95 ± 2%) whereas the poloidal 
field tended to be more non-axisyme tric (65 ± 7%). This 
is con siste nt with the observation s of iDonati fc Landstreetl 
l|2009l) and IMarsden et al] l|2011ah where stars with signifi- 
cant toroidal fields often have non-axisymetic poloidal fields. 
The error estimates on these measures were obtained by 
varying a range of stellar parameters such as inclination, ra- 
dial velocity, f2eq, 5^, period andiisini. Hence these may be 
refered to as a variation bar, which is analogous to an er- 
ror bar except that the measurement indicates the variation 
found by this method. 

When measuring the differential rotation using the mag- 
netic features, the equatorial rotation rate iloq was 4.51 ± 
0.01 radd~^, with a photospheric shear 50, of 0.24 ± 0.03 
radd-^ The comparison with those parameters found using 
the Stokes / information is shown in Figure 1121 However, 
this difference in the 50 values are within the error bars 
of the two measurements. When compared with the differ- 
ences observed on other stars of simil ar spectral type such 
as HD 141943 (|Marsden et al.ll2011bl) . these features may 
be anchored at similar depths in the convective zone, unlike 
HD 141943 which may be anchored at different depths of the 
convective zone. However, this conclusion must be treated 



with caution until further Stokes V data is obtained for this 
star. 



4.3 Hq emission 

Figure [9] shows a bright spot in this dynamic profile at ~0.15 
phase. This is where the chromosphere is much denser as a 
result of an increase in temperature in that local region. 
There appears to be a reduced emission at ~0.75 phase as 
well. These variations in Ha emission lie within the stellar 
velocity profile and hence arise from regions close to the stel- 
lar surface. We thus interpret these Ha variations in terms of 
prominences more akin to solar prominences than the more 
extended features knows to exist at the co-rota tion radius of 
other active solar-type stars, such as LQ Lup (|Donati et al.l 
|2000D . 



5 CONCLUSIONS 

In this paper, reconstructed brightness and magnetic images 
of the young solar-type star HD 106506 have been presented. 
The brightness image shows low-latitude to mid-latitude fea- 
tures coupled with a predominate polar spot. This is very 
similiar to many of rapidly rotating solar-type stars. The 
magnetic images reveal regions of azimuthal field near the 
stellar surface that suggest that the dynamo mechanism may 
be occuring entirely in the convective zone and close to the 
surface of the star itself. The p hotospheric shear f or HD 
106506 supports the findings of iBarnes et all |200i), that 
effective temperature is the dominant factor in the level of 
surface differential rotation. However, the reason why such 
a relation exists remains unclear. 
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Table Al. Transformation Coefficients for the 25th July, 2005 
photometry at MKO. The associated errors, rms fit and value 
for each filter have been included 



n 











-3.217 


-3.506 


-3.709 




±0.14 


±0.13 


±0.12 


1 


0.311 


0.144 


0.134 




±0.11 


±0.09 


±0.10 


2 


-0.154 


0.065 


0.020 




±0.12 


±0.06 


±0.14 




0.103 


0.060 


0.023 


rms fit 


0.021 


0.019 


0.011 



APPENDIX A: BVR PHOTOMETRY AT MKO: 
2005 AND 2007 

The IRAF package was used for the reduction and analy- 
sis of the data. The CCD frames had the dark current re- 
moved (that also contained the Bias information) and were 
flat fielded using a nightly master flat field, which was a com- 
bination of 9 dome flat fields. To produce a single data point, 
three images taken through each filter of each star field were 
averaged, using IMCOMBINE, to produce an image with 
improved signal-to-noise from which the magnitude of all 
the program stars were measured. The stars in the frame 
were found using DAOFIND and simple aperture photom- 
etry was done using DAOPHOT. The sky background was 
determined and removed using the mode of an annulus of 10 
pixel radius a round the star. The curves of growth method 
l|StetsoDlll990l ) was used to determine the aperture size that 
would be used to determine the magnitude of each of the 
stars in the field. 




Figure Al. The area surrounding HD 106506, showing 5 of the 
brighter stars were used in the determination of the period. The 
stars' names are shown in Table lA2l in appendix A. 



Table A2. The BVRI information for HD 106506 and the com- 
parison stars. Graham (1982) E-Regions, E5,E6 and E7, were used 
to standardise these stars and the data was processed using the 
routines in PHOTCAL found in IRAF. The data were taken on 
the 25th July, 2005 and a sample image, with the star number 
included, is shown in Figure IAT] 



Star 


Star's Name 


V 


BV 


VR 


Target 


HD 106506 


8.463 


0.594 


0.356 


Star 1 


HD 106507 


9.152 


-0.024 


0.044 


Star 2 


TYC 8979-979-1 


10.108 


0.100 


0.101 


Star 3 


TYC 8978-1180-1 


10.003 


1.083 


0.584 


Star 4 


HD 106417 


10.133 


1.015 


0.558 


Star 5 


TYC 8978-5479-1 


10.948 


0.203 


0.165 



Al Standardisation of Comparison Stars: 25th 
July, 2005 

BVR photometry was conducted at MKO in 2005 in or- 
der to determine the photometric colour of HD 106506 and 
the other nearby stars within the field of view. Figure lAll 
shows the target star and 5 comparison stars within the 
same field of view, are indicated. Several BVR photometric 
observatio ns were taken o f the standard stars in Graham's 
E-regions (ICrahaml Il982l ). The O'Mara telescope used for 
this procedure comprises a 35cm Schmidt-Cassegrain optics, 
an SBIG STL-1301E CCD camera with 1280x1014 16 
square pixels, and a Paramount ME mount. The three re- 
gions were used, the E-5 (Stars A,0,U,S,V,Y,c), E-6 (Stars 
98,M.P,X,W) and E-7 (Stars M and S) regions. The rou- 
tines found in the PHOTCAL package in IRAF were used 
to determine both the transformation coefficients and the 
magnitude of the stars, including HD 106506. There were 
three images taken through each filter of the star fields and 
then averaged using IMCOMBINE, to produce a single im- 
age from which the magnitude of all the standard stars and 
the program stars were measured. The following transfor- 
mation equations were solved using the FITPARAM task in 
IRAF: 

mb = B + 6o + 6iAb + b2(B-y) (Al) 



mv = V + ?;o + i'iAv + t)2(B- V) (A2) 

nir^ R + ro + riX, + r2{V - R) (A3) 

where B, V and R are listed magnitudes of the stan- 
dard stars in lGrahaml(|l982l ). The instrumental magnitudes, 
as measured by DAOPHOT were mb, rriv and rrir. The zero- 
points are bo, wo and ro; while b\,vi and ri are the extinction 
coefficients and 62,^^2 and r2 are the colour correction coef- 
ficients. The values determined for the coefficients are given 
in Table lAll and the results are shown in Table IA2I 

A2 BVR Broadband Photometric Observations 
from MKO: 2007 

Table[T]shows a log of data taken during the observation run 
in 2007. The aperture size selected had a pixel radius of 15, 
giving an aperture radius of approximately 8 arc seconds. 
On the 11th of April, several observation s were taken of the 
standard stars in Graham's E5-regions ()Graham|[l983 ) so 
as to place HD 106506 on the standard magnitude scale. 
The stars used were A,0,U,S. The routines found in the 
PHOTCAL package in IRAF were used to determine both 
the transformation coefficients and the magnitude of the 
stars, including HD 106506. There were three images taken 
through each filter of each star field and these were aver- 
aged, using IMCOMBINE, to produce an image with im- 



Magnetic fields and differential rotation on the pre-main sequence. 
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Figure A2. The top left panel shows the V-band light 
curve for HD 106506, with the phase referenced to epoch = 
2453543. 93642361-d which coincided approximately with the mid- 
dle of the AAT run. The BV and VR colour light curves are also 
shown. The bottom right panel shows the Discrete Fourier trans- 
form, which is interpreted in terms of a rotation period of 1.416 
days. 



proved signal-to-noise from which the magnitude of all the 
standard stars and the program stars were measured. 



A3 Period Determination in 2007 

The rotational period was determined using Discrete Fourier 
transform s (DF T), using the algorithm as explained by 
iBelseren^ (|l988l l. The five brighter stars were selected as 
possible comparison stars and the period was determined 
for each in B, V and R band. The average period was deter- 
mined to be 1.416 ± 0.133 days, and was consistent with all 
the comparison stars in B, V and R bands. The error was de- 
termined by finding the full width half maximum of the DFT 
line profile. Figure IX2l shows the phased folded light curve 
for the V- band as well as the BV and VR colour indices. 
The power spectrum is also shown as generated by the DFT 
algorithm. Graham's E5 Region, coupled with the 5 stars 
previously standardised, were used to place HD 106506 onto 
the Cousins magnitude system. 

This paper has been typeset from a TJjX/ M^jX file prepared 
by the author. 



